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Abstract 

The dynamics of a single hole in quantum antiferromagnets is influenced 
by magnetic fluctuations. In the present work we consider two situations. 
The first one corresponds to a single hole in the two leg t — J spin ladder. 
In this case the wave function renormalization is relatively small and the 
quasiparticle residue of the S = 1/2 state remains close to unity. However at 
large t/J there are higher spin (S = 3/2, 5/2, ..) bound states of the hole with 
the magnetic excitations, and therefore there is a crossover from quasiparticles 
with S = 1/2 to quasiparticles with higher spin. 

The second situation corresponds to a single hole in two coupled anti- 
ferromagnetic planes very close to the point of antiferromagnetic instability. 
In this case the hole wave function renormalization is very strong and the 
quasiparticle residue vanishes at the point of instability. 

PACS: 71.27.+a, 75.10.Jm, 75.50.Ee 
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I. INTRODUCTION 



Quantum spin ladders have recently become a subject of considerable interest, mainly 
because of their relevance to a variety of quasi one-dimensional materialsS. In addition 
they represent, on the theoretical side, an important step towards understanding of the 
two-dimensional Cu02 planes in the cuprate superconductors. The properties of a single 
hole in the two leg ladder were investigated numerically by Troyer, Tsunetsugu, and Rice! 
and Haas and DagottoQ using exact diagonalizations. More recently Edera approached the 
same problem using both exact diagonalization and the Self Consistent Born Approximation. 
Two most interesting properties have been observed in these works: The first one is the shift 
of the hole dispersion minimum to the point close to k « n/2. The second property is the 
closeness of the quasiparticle residue to unity even at large t/J. Such behavior of the residue 
is very much different from the t — J model on the 2D square lattice!. 

The closeness of the quasiparticle residue to unity indicates that a simple variational 
approach and even perturbation theory must be applicable to the problem, and therefore we 
use these techniques in the present work. This allows us to find the 5 = 1/2 quasiparticle 
dispersion analytically. The dispersion is in an agreement with previous computations. We 
demonstrate also that there are higher spin (S = 3/2, 5/2, ...) bound states of the hole with 
the magnetic excitations. At large t/J these bound states lie below the S = 1/2 state, and 
this means that they become the real quasiparticles of the system. So our conclusion is 
that the spectrum of the single hole states on the ladder is very rich. This complexity is 
somewhat similar to the complexity of the spin-wave spectrumEm. 

In spite of the complex spectrum, the structure of each individual state on the ladder is 
relatively simple. It means that the number of significant components in the wave function 
remains of the order of unity. The purpose of the second part of this work is to consider a 
spin liquid with really complex structure of the hole wave function. We demonstrate that 
this situation is realized in a 2D spin liquid close to the point of antiferromagnetic instability. 
Specifically we consider two coupled antiferromagnetic planes. As the point of instability 
is approached the number of components in the hole wave function approaches infinity and 
the quasiparticle residue vanishes. 

II. T - J LADDER 

We consider the standard t — J model on a two-leg ladder. The Hamiltonian reads 

H = H t + Hj, (1) 
H t = - ]T (t c\ a c Ja + H.c) , 



<i,i> 

where cf CT is the creation operator of an electron with spin a at site i. As usual double 
occupancy is forbidden, (ij) represents nearest neighbor sites. The spin operator is Sj = 
^cl a a a pCip. For bonds along the legs we choose = J, and for bonds along the rungs 
Jij = J±- We are mostly interested in the "isotropic" case: J = J±, but it is convenient 
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to keep J± > J as an independent parameter. Following Chubukovi, Sachdev and Bhatti 
and Gopalan, Rice, and Sig rist0 we consider the ground state consisting of inter-chain spin 
singlets \GS) = |1, s)\2, s)|3, s)..., where \n,s) = ^= [| |)i - | T)j] as a starting point 
of our approach. Since each singlet can be excited into a triplet state it is natural to introduce 
a creation operator tf an for this excitation: 



K«) = 4nK s )> a = x,y,z. 



(2) 



The interaction of the triplets between themselves is described by the Hamiltonian (see, e.g. 
RefsH) 



J 



+ 7 («U 



f3n+lt/3ntan+l 1 an L an+l L frn L (3n+l 



In addition, the hard core constraint 







(3) 



(4) 



has to be enforced on every sited. The interactions (§,f|) lead to rearrangement of the 
unperturbed ground state and nontrivial spectrum of the spin-wave excitations0'i'i. 

If we remove one electron from rung n a single hole state on this rung is created. Let 
us denote by a\ ima the creation operator of a hole with spin a on rung n, where the index 
m — 1, 2 represents the legs of the ladder. This operator is defined by the equation 



a ^n,m,a\ n i S ) ~ C n,rn,a\fy i 



(5) 



where m = 2 if m = 1, and fh = 1 if m = 2. This is a fermionic operator which satisfies 
the usual fermionic anticommutation relations. The hole on the rung n interacts with the 
triplet excitations on its nearest sites. This interaction can be easily found by calculating 
all possible matrix elements of the initial Hamiltonian (|l|) (see Ref.S). The result is 



H h = ~tJ2 ( a n,2,a W + H - c ) + « £ (' 



a n+l,m,cr a n,m,a 



+ H.c. 



(6) 



+ e k 



t 



t Sn±i,n + 2 S n= i-i )n -i-i 



n,m,a 

H.c. 



+ r, E (^n+l^n a i,m,o- a n+l,m,o- + -^' C - 



n,m,cr 

2 S n +l,n 



X t 



n+1 



+ ff.c.) + -r E S "' n t ™ ±1 X tn±1 + H ' c 



n,± 



The index ± denotes summation over nearest rungs. Following Ref.@ we have introduced 
the notations 



m,a,/3 



(7) 
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where a is the vector of Pauli matrices. In addition we have to impose the constraint that 
a triplet excitation and a hole can not coexist on the same rung 

The effective Hamiltonian given by Eqs.(|3]) and (H) (H e ff = H sw + H h ) combined with 
the constraints (|4]) and (|) represents an exact mapping of the initial Hamiltonian ([!]) in the 
sector with a single hole". One can work with H e fj by using standard diagram techniques 
in momentum representation. The constraint ([8]) can be taken into account in a mean field 
fashioni or by introducing a fictitious (ghost) interaction similarly to the way it was done 
for spin waves in the paper0. However fortunately the problem can be resolved in a simpler 
way. Because of the sizable spin-wave gap there is no real long-range dynamics on the ladder. 
Therefore the problem can be solved in the coordinate representation by constructing a trial 
wave function that explicitly obeys the constraints (|j) and (H). First consider an isolated 
rung with a hole. Say this is the rung n. The energy of this state with respect to the 
undoped system is 

3 3 J 2 

* = 5 Jx + - iT± . (9) 

The origin of the first term in eo is obvious: one antiferromagnetic link is destroyed. The 
origin of the second term is also simple. There are quantum fluctuations in the spin liquid 
state: virtual excitations of the triplet pairs created by the term ^t^+i in the Hamiltonian 
(H). This gives vacuum energy 3 P er nn ^' wnere the coefficient 3 is due to the number 
of possible polarizations!!!. The presence of a hole on rung n suppresses fluctuations at the 
links n—l,n and n, n + 1 and this gives the second term in ([]). Because of the hybridization 
along the rung the stationary states are symmetric and antisymmetric with hopping energies 
—t and +t respectively. Thus for an isolated rung one finds 

~ a "' 1 ' <J + a "' 2 '- , e = eo-t, (10) 



f _ "71,1,(7 a n,2,a _ , 

a n-,a — ^J= j C — Co + t. 

Consider now a hole hopping from rung n to n + 1: aJ n + ♦ — > a^ +1 + *. Naively, according 
to the Hamiltonian @ the effective matrix element for this process is t/2. However this 
does not take into account quantum fluctuations of the spin-wave vacuum. By taking into 
account these fluctuations we obtain the following initial wave function 
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, 3 j- + 

v l — 2 A* ) + ^a,n-2^a,n-l 



at 



(1 — 3/i 2 ) + ^^,n+l^,n+l + A t * J 9,n+2* l 9,n+3 I s )' (H) 



Here we explicitly present the part of the wave function corresponding to the rungs from 
n — 2 to n + 3. The coefficient \i describes the spin- wave quantum fluctuations and according 
to the Hamiltonian (|3|) it is equal to 

„— ^ (12) 
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The final state is 



I/) — (1 3/i 2 ) + [lt a n _2^L n-l + A^cc, n-l^t n 



r/ 



t 

n+l,+,t 



, 2\ t t 

J- — 7^ U ) + A^/3,n+2^/3,n 



+3 



Calculating the matrix element of the Hamiltonian H t , Eq. (|T]) we find 



:P , P 



2' ' ' 2' 

which agrees with that obtained in the Reffl 

Let us consider now the possibility of a virtual spin-wave excitation 



\s). (13) 



(14) 



V+.T 



f f 4 



(15) 



which gives a second order correction to the hole energy. Note that the symmetry of the 
operators a) in the admixed components is negative. The Hamiltonian ([!]) is symmetric under 
permutation of the ladder legs. Therefore all excitations can be classified according to this 
symmetry. The operator v has (— )-symmetry, and therefore the a) which appears together 
with v must be of the (— )-symmetry as well. Below in the notations of the quantum states 
we will omit the symmetry index. The energy difference corresponding to the excitation 
©is 



AE 



e + t + J ± 



1 3J^ 
3 + 8J, 



(e -t) = 2t + J± 



3J! 
8J 



(16) 



The admixture (JT5|) first of all shifts the on site energy. Using the wave function ([TTj) we 
find 



(4 a n+il#tl4> " „ (4+i4I#jI4> 
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3pH 2 3p 2 J 2 



AE AE 2 AE 

The admixture (|T3j ) generates also an effective one step hopping 



8 AE 



(al +1 \Hj\ ^4+i ) (t n 4+i |# t |4) (4+i 1 H t \ t\ l+l al) (£ n+ i4 | #j 1 4) 



At 



ft 



AE 



(17) 



(18) 



3 p 2 tJ 
4 ~AE' 



There is one more contribution to the effective one step hopping which arises due to the spin- 
wave quantum fluctuations: Let the hole first reside on rung n, then a spin-wave quantum 
fluctuation arises on rungs n+1 and n + 2, after that the hole hops to rung n+1 annihilating 
the spin-wave excitation t n+1 and after that the spin wave excitation t n+2 is annihilated 



because of the interaction Hj with the hole a 



t 

n+1- 



There also exists the time reversed 



mechanism which gives an equal contribution. This gives 
-2fi 



(4+1 \Hj\tn+2 a n+l) (4i+2 a n+l 



-^t|^n+2^n+l4) 



AE 



3 \iptJ 
I AE ' 



(19) 



A similar mechanism generates a two step hopping: Initially the hole is on rung n, then 
a spin-wave quantum fluctuation arises on rungs n + 1 and n + 2, and after that the hole 
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annihilates this virtual excitation moving to the rung n + 2. Again the time reversion of this 
mechanism leads to an equal contribution. 

,(2) _ _2 ( a n+2 1 Ht I ^n+2 a n+l ) (^L+2 Q n+l I ffit I ^n+2^n+l fl n ) _ 3 /ipt 2 

Adding together (|T7|) , (pT8[) , (|I9"D , (f20|) we obtain the total second order correction to the hole 
energy 

Sef = [pH 2 + ip 2 J 2 + t J(p 2 + Pf j) cosg + 2ppt 2 cos2g) . (21) 

This formula, combined with ( |T0| ) and ( H ) gives a reasonable description of the hole disper- 
sion because the energy interval AE is large enough. 

The accuracy of the calculation can be improved by replacing the perturbation theory by 
a variational method. The zeroth approximation hole wave function and the corresponding 
dispersion are of the form 

^t = ^E^ n < +lT > (22) 



ei = eo — t + p 2 t cos q. 

The admixed components represent the hole and the spin wave combined into total angular 
momentum | with z-projection | 

l+) = ^E^ n (4<4fi)„. 

|-> = ^£^(*W-i) ai > (23) 



'N 

e 2 = ei + AE. 



The mixing matrix elements are 



-2i<j 



V = (+|i?|^ T > = h) = " ^P Je ' iq ~ \^ Q - (24) 

Diagonalization of the matrix gives the quasiparticle dispersion and residue (weight of the 
component (|22|) in the total wave function) 



e 



e l + AE/2-^/(AE/2Y + 2\V\ z (25) 



Z=i|l 



2\ </l + 8|Vf /(A^) 2 



Expansion of this solution in powers of (V/AE) gives the perturbation theory result (21). 

We are interested in the "isotropic" ladder, Jj_ = J. In this case fromperturbation 
theory ([12]) , p = —0.25. Comparison with more accurate computationsoao shows that 
really (i = —0.3. The difference is small, but nevertheless in further numerical estimations 
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we will use \x = —0.3 which effectively accounts for some higher order corrections. At t <C J 
Eq.(p5|) gives 

e q « 1.28 J- 0.7t + 0.24tcosq, (26) 
Z = 0.85 

The plots of the dispersion for tj J — 1, 2, 3 are given in Figures 1,2,3. They agree with results 
of numerical computationsi~i. The quasiparticle residue is close to unity. For example for 
t/J = 2 the residue Z « 0.9 at q = 0, Z « 0.83 at q = vr/2, and Z « 0.97 at g = vr. 
To avoid misunderstanding we note that this residue can not be directly compared with 
"photoemission spectra" because of two reasons. First, the "photoemission" operator is 
different from (see the discussion in the Ref B) . The second reason is that the quasiparticle 
we consider is well defined only near bottom of the dispersion. At small q the quasiparticle 
energy is large and therefore it is unstable with respect to real emission of a spin wave which 
we neglect in our consideration. 



III. HIGHER SPIN HOLE STATES IN THE LADDER 

It is clear that at very large t the hole tends to order all the spins ferromagnetically. 
This is the Nagaoka effect, as discussed in the Ref.!. The same effect produces higher spin 
bound states of the hole with the spin wave. Let us consider first the state with S = 3/2. 
The simplest ansatz for this state is 



where t^ +1 = ;4g(— 4, — it\) corresponds to spin projection +1. The energy corresponding to 
this wave function is 

e^=e -2t + J ± + i + l^- (28) 

The first term is the single hole magnetic energy the second term is due to hybridization, 
J± is the excitation energy of the triplet, J/4 is due to the antiferromagnetic interaction 
between the hole and the triplet, and the last term is due to the blocking of quantum 
fluctuations. The state ( |2"7| ) can propagate via intermediate virtual configurations: first the 
triplet hops and then the hole catches up. This gives an effective hopping matrix element 

,(3/2) tJ / A (29 ) 

W/ " t + 3.77(8.70 - J/4" [Zy) 

Thus the dispersion of the spin 3/2 state is 



e 



(3/2) _ (3/2) , 9 .(3/2) 



<-l 



W + 2t%;><x*q. (30) 



The plots of (p0|) for tj J = 1, 2, 3 are given in Figures 1,2,3. 
Finally, let us consider a ferromagnetic bag with S = 5/2 
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^ — /g (^+l,n^+l,n+l a n+2,+,t + ^+l,n a n+l,+,t^+l,n+2 + ^+l,n^+l,n+l a n+2,+,y ■ (^1) 

The energy of this state (neglecting the dispersion) is given by 

eOW = eo - (1 + >/2)t + 2J ± + ^ + (32) 

4 4 J_|_ 

For t/ J = 2, 3 it is plotted in Figures 2 and 3. 

Now let us discuss our results. In the accepted approximation the bag state with spin 
S = 3/2 at arbitrary t/J has higher energy than the bottom of the quasiparticle (S = 1/2) 
dispersion. Therefore the S = 3/2 state is unstable with respect to decay into a quasiparticle 
and a spin wave". It exists only as a resonance in the quasiparticle - spin wave scattering. 
The 5=1/2 quasiparticle is itself unstable in the vicinity of q — 0. It can decay into a 
quasiparticle and a spin wave with higher q, or into an S = 3/2 state and a spin wave. 

The energy of the S = 5/2 bound state decreases with t faster than that of the S = 1/2 
state. Therefore sooner or later the S = 5/2 bag state will become the lowest one and thus 
the true quasiparticle of the system. Comparing (p5|) with (|32|) we find the crossover point 
t c « 7 J. As t increases further, in accordance with Nagaoka's theorem, quasiparticles with 
higher spins will appear. 



IV. SINGLE HOLE IN TWO COUPLED ANTIFERROMAGNETIC PLANES 

The Hamiltonian in this case is similar to ([I]). The only difference is that instead of 
the two-leg ladder we consider now two coupled planes (we assume that both planes form 
square lattices). The coupling J is the antiferromagnetic exchange in the plane, and J± is 
the antiferromagnetic exchange between the planes. It is known that at J± > J c ( J c ~ 2.5 J) 
the system is in a disordered dimerized spin liquid state, and at J± < J c it undergoes 
quantum phase transition to a Neel phase with long range antiferromagnetic order (see, e.g. 
Refs.00). In the present section we will consider the spin liquid phase in the regime very 
close to the point of antiferromagnetic instability. 

Let us remove one spin from the system. We will see that the most interesting situa- 
tion which is qualitatively different from the ladder arises when the hole is static (t — 0). 
Therefore, let us consider the static hole (impurity) problem. In the first approximation, 
neglecting the interaction with spin waves, the energy of the hole is 

eo = 5jj. + ~«2.5J. (33) 

It is similar to @. The first term is due to the destroyed bond between the planes, and the 
second one arises from the suppression of quantum fluctuations on the four links in the plane. 
We substitute in ([33D J± = J c ~ 2. 5 J. Interactions with spin waves give corrections to (j33|). 
The hole-spin- wave vertex is represented by the vS term in the effective Hamiltonian (^j). 
In momentum representation it has the form 

9k,d = -^7q(Mq + Wq)<X a (34) 



S 



This vertex describes the process when a hole with momentum k decays into a hole with 
momentum k — q and a spin wave with momentum q and polarization a, (a = x, y, z). The 
spin Pauli matrix of the hole is denoted by a; 7 q = \{cosq x + cosg^); w q and t> q are the 
Bogoliubov parameters corresponding to the spin wavetil. In ([34]) we neglect the deviation 
of the spin wave quasiparticle residue from unity because it gives only a 10% correction 
sec Ref.0). The spin wave dispersion has a minimum at Q = Near this point the 

dispersion is 



Wq = ^A2 + c 2( q -Q)2, (35) 

where c ~ 1.9J00 is the spin wave velocity. At q ~ Q the Bogoliubov parameters can be 
represented asQ 



UqWVqWi/ , (36) 



v 2w q 

where A ps 2.4J. 

The naive second order correction to the hole energy is given by 

Se = -}^ a = — —5- / — « - J In — , 37 

where the integral is calculated with logarithmic accuracy. We have introduced here the 
notation 

3J 2 A 

J = =- « 0.63J (38) 

7TCT 

It is clear that since ^ 3> 1, the naive formula ( j3"7|) has no meaning. However, the presence 
of the logarithmic divergence indicates that the self consistent Born approximation (see e.g. 
Ref.Ej'L<J) i s valid in this case. Following this approximation let us consider the quasiparticle 
(pole) part of the hole Green's function 

G{e) = Z — { -, (39) 

f-e+tO 



where e is the energy of the state. This Green's function must obey Dyson's equation 

1 

eo - S(e) 

where the self energy is given by 



G{e) = (40) 



m = 3 E <q G(e - Wq ) « Z / — - *± — ■ (41) 



2tt 2 J u; q (e — e — u; q; 



The calculation of this integral is straightforward and gives 



Z(e) = -ZJ ln± f^E(e)) = (42) 
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Finally, substituting into ((E]) and comparing with (p9|) we find 

e = e - ^AJo In ^ (43) 




The quasiparticle residue Z vanishes at A — > 0. It is clear that it does not happen at finite 
hopping, because in this case the lower limit of integration in ( |37|) is t and the integral 
remains finite at A = 0. Physically it means that due to the finite speed of propagation the 
hole has no time for dressing by an infinite number of spin waves. 

Thus, the static impurity is strongly dressed by spin waves. The question arises: is this 
cloud dilute enough to justify the spin-wave approach we have applied, or is it dense in 
which case the above consideration has no quantitative meaning because we do not take 
into account the interaction of the spin waves between themselves? The size of the cloud 
is r ~ c/A (in units of the lattice spacing), and the effective number of spin waves is 
N sw ~ 1/Z — J Jo/ A. Therefore the density of spin waves n ~ N sw /r 2 ~ Y_^o^3/2 rema i n s 
low and hence the dilute gas approximation is valid. 

V. CONCLUSIONS 

Using variational approach and perturbation theory we have found analytically the dis- 
persion of a single hole with spin S = 1/2 in the t — J ladder. The dispersion is in an 
agreement with previous computations. We demonstrate also that there are higher spin 
(S = 3/2, 5/2, ...) bound states of the hole with the magnetic excitations. At large t/ J these 
bound states lie below the S = 1/2 state meaning that they become the real quasiparticles 
of the system. We have estimated the critical value of t/ J where the crossover to the higher 
spin quasiparticles take place. 

In spite of the complex spectrum, the structure of every state on the ladder is relatively 
simple. This means that the number of significant components in the wave function remains 
of the order of unity. In the second part of the present work we have considered the opposite 
situation when the wave function has a very large number of significant components. This is 
realized for a static hole in two antiferromagnetic planes close to the point of antiferromag- 
netic instability. We demonstrate that the self consistent Born approximation is justified 
in this case. Using this method we show that as the point of instability is approached the 
number of components in the hole wave function approaches infinity and the quasiparticle 
residue vanishes. 
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FIGURES 



FIG. 1. Dispersions of the 5=1/2 quasiparticle (solid line) and the 5 = 3/2 bound state 
(dashed line); t/J = 1. 

FIG. 2. Dispersions of the 5 = 1/2 quasiparticle (solid line), 5 = 3/2 bound state (dashed 
line), and 5 = 5/2 bound state (dashed-dotted line); t/J = 2. 

FIG. 3. Dispersions of the 5 = 1/2 state (solid line), 5 = 3/2 state (dashed line), and 5 = 5/2 
bound state (dashed-dotted line); t/J = 3. 
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